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Abstract. We study the asymptotic behavior, as Planck's constant H — > 0, 
of the number of discrete eigenvalues and the Riesz means of Pauli and Dirac 
operators with a magnetic field /i~B(x) and an electric field. The magnetic field 
strength /i is allowed to tend to infinity as h — * 0. Two main types of results 
are established: in the first /ih < constant as h — * 0, with magnetic fields of 
arbitrary direction; the second results are uniform with respect to [i > but 
the magnetic fields have constant direction. The results on the Pauli operator 
complement recent work of Sobolev. 



1. Introduction 

The Dirac and Pauli operators Dy,P^i which are the objects of study in this 
paper, are defined as follows : 

Dy = Dy(B) := a ■ (tV - fM) + + V 

i 

3 

(1) = ^a k {-ihd k - LM, k )+ p + V 

fe=i 

P H /eeP h/ (B) := [a ■ (-ihV - fia)} 2 + W 

3 



[^2a k {-ihd k -fia k )} 2 + W 



k=l 

(2) = H (B) - fiTur ■ B + W, 

where Hq(B) — (— ihV — /ia) 2 is the Schrodinger operator with magnetic field B 
and 

• a = (ai, d2, as) is a magnetic vector potential with magnetic field B := V x a; 



a = (<7i, ct 2 , (T3) is the triple of Pauli matrices 

i\ fo -A (\ 



r>1 = \\ o)> a2= \i 0)' CT3= vo -1 

a = a.%, a-i) where oej, j = 1, 2, 3, and (3 are the Dirac matrices 
( O2 <7k\ , , r, o a (h 2 

afe= U o 2 J' fc = 1 ' 2 ' 3 ' P={o 2 -I 2 

where O2 , 12 are the 2x2 zero and unit matrices respectively. 
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Given a domain ft C R 3 , self-adjoint realizations Piy(ft) and Py(O) are defined 
by Dirichlet boundary conditions on ft, these being satisfied in the usual weak sense 
(see §2). The Pauli operator ¥ W (Q) acts in L 2 (ft) ® C 2 = [L 2 (tt)] 2 and the Dirac 
operator By (ft) in L 2 (ft) ® C 4 = [L 2 (ft)] 4 . Under quite general conditions the 
spectrum of Po coincides with [0, oo), and the perturbation W introduces negative 
eigenvalues A„(Pw(ft)). The Dirac operator D has typically a spectrum R\ (— 1, 1) 
and V causes eigenvalues A„(IV(^)) to appear in the gap (—1, 1). Our concern in 
this paper is with the Riesz means 

(3) M 7 (P w ,ft) = ^|A„(P w (ft))p\ M 7 (Dv,n) = 5^|An(10v(n))r 

n n 

where 7 > 0, and, in particular, the counting functions given by 

(4) N(P W , ft) = M (P W , fi), N(B v ,n) = M (B V , ft). 
In [p3f , Sobolev investigated the natural quasi-classical formula 

(5) M 7 (P W , ft) - ft~ 3 <8 7 (^|B|, W, ft), ft -> 0, 
with the "magnetic" Weyl coefficient 

oc 



(6) Q3 7 (6,w,ft) := /3 7 / 6(x) 



fc=i 



where 



/ 3 7 = i ^^ 1 ^(i-r 1/2 *. 

Under weak regularity conditions on B and W, Sobolev established two main re- 
sults: 

(i?i): for fih < constant, (||) is satisfied for 7 > 1; 

(-R2): if B has constant direction, (||) is satisfied uniformly in ^ > for 7 > 1/2. 
Earlier results of this type were proved for homogeneous (i.e. constant) fields by 



Lieb, Solovej, and Yngvason 15 , |16[, and for non- homogeneous fields by Erdos 
and Solovej (tJ, ||, with /ih 3 —* 0. Our objective in this paper is to investigate the 
validity of (i?i) and (R2) for all 7 > 0. We prove in Theorem [l] that if W, |B| S 
L 3 / 2 (ft) in (i?i), then (||) is satisfied for all 7 > 0. This extends a result in || in 
which a Weyl asymptotic formula is established for the case fih ~ as K — * 0. If 
|B| does not belong to L 3 / 2 (ft) in (i?i), we obtain the result for N(P W + A, ft), the 
number of eigenvalues of P^(ft) less than — A < 0. Our main result in problem (R2) 
is also of this form. This is the best that can be expected in general for 7 = in 
view of the absence of an inequality of Cwikel, Lieb, Rozenblum type for the number 
of negative eigenvalues; such an inequality is available when |B| S L 3 / 2 (ft). We 
also investigate circumstances in which (||) holds for 7 > when B has a constant 
direction. 

In jl8) the leading term in the asymptotic value of N(Pw + A, M 3 ) is determined 
for A > 0, H= 1,B of constant direction, and /1 — > 00. Specifically, with B(x) = 
(0, 0, B(x)), x = (x, X3), x G K 2 , where B(x) is bounded above and away from zero 
and has a bounded gradient, 



(7) lim fi^NiPw + A,R 3 ) = D(A) := ^- / N(X 



w+\(x),R)B(x)dx, 
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where, for any fixed x £ R 2 , Xw(x) is a self-adjoint realization of —4^ + W(x, •) 

ax 3 

in L 2 (R) with essential spectrum [0,oo). From many other results on the spectral 
asymptotics of N(¥\y,Cl), we mention in particular that of Iwatsuka and Tamura 
in 01 for h = n = 1 and A -> : 

(8) N(P W + \,M. 3 ) ~V(\). 

Results for the two-dimensional problems are also obtained by the authors cited 
above, and our techniques also apply in this case. Furthermore, we derive analogous 
results for the Dirac operator with magnetic field in both cases (R%) and (R2) above. 

Our strategy is based on that in p3| which in turn was inspired by ideas from 
§). For the case fih <constant, the technique we use involves a tesselation of R 
by cubes Q, the derivation of two-sided estimates for N(P\y,Q), with constant 
W and B, in terms of the explicit values of the eigenvalues (Landau levels) of 
the operator realization of ¥w on & torus, and subsequently the localization of 
V w (^l) in terms of the Pw{Q) an d the application of either the Cwikel, Lieb, 
Rozenblum (CLR) inequality for the magnetic Schrodinger operator (in the case 
in which there is known to be a finite number of negative eigenvalues in problem 
(i?i)), or an inequality derived from a Lieb-Thirring inequality established by Shen 
in p(J for the trace of the negative eigenvalues. For problem (R2) this method is 
too crude. Sobolev in ]23] used the spectral properties of Fyy on a torus with an 
arbitrary periodic magnetic field B having integer flux, and, in particular, that zero 
is an eigenvalue of the operator with multiplicity equal to the flux of B. Another 
important ingredient in |23| in this case is a Lieb-Thirring inequality established in 
(22J in which |/zB| has the proper (linear) scale; for a discussion of the significance 
of this see the paper of Erdos and Solovej @, §1. We use a similar estimate 
derived by Shen in |20[] for N(¥\v + A,R 3 ) with A > 0. This grows like 1/y/X as 
A — » 0, and implies (||) only when 7 > 1/2, as already established in [^3| with 
similar assumptions. We also present a result for a cylindrical domain Q in which 
N(¥w + A, f2) grows logarithmically, which implies that M 7 (Pw, is finite for all 
7 > (see Proposition ||). 

Throughout the paper (•, -)n, || • ||n are used to denote the usual inner-product 
and norm in each of the Hilbert spaces L 2 (f2), [L 2 (f2)] 2 and [£ 2 (f2)] 4 ; the precise 
space will be clear from the context. We adopt the convention that inner-products 
are linear in the second argument and conjugate linear in the first. We denote the 
identity matrix on C 4 = C 2 C 2 by I 2 . We denote points in R 3 as x = (x, x 3 ), x £ 
R 2 , X3 £ R. We shall write A<B, or B>A, to mean that \A\ < CB for some positive 

constant C, and /(x) = 0(A) will mean that /(x) is bounded by A. 

2. Preliminaries 

To begin, we need to define the operators precisely and locate their essential 
spectra. The Pauli operator will be defined as the form sum of P and the 
operator of multiplication by W, the conditions assumed on W being sufficient for 
the associated sesquilinear form to be lower semi-bounded. For the Dirac operator 
we require the restriction of By to Cq°(&V) to be essentially self-adjoint, in which 
case EV is the self-adjoint operator associated with the form ||By0||n. In §3 below, 
our requirements in Theorem |l| are met by the following results from Define, 
for m = 2,4, 

Hi = nl{R 3 ) := {u : u, {{h/i)d„ - fia v )u £ [£ 2 (R 3 )]"\ v = 1, 2, 3, } 
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in which the derivatives are defined in the distributional sense. If a S [£ 2 oc ] 3 , 
then [C£°(R 3 )] m is dense in H\ (see Kato Q, Simon @, and Leinfelder and 
Simader Mj). Also, the diamagnetic inequality 

|(fi/»)^(|tt|)| < \{{h/i)d v -iia v )ul u= 1,2,3, 

holds for almost every x £ R 3 and u G (see Lieb and Loss jl3), p. 179). As a 
consequence, u i— ► |u| maps Ti} a continuously into the Sobolev space iJ 1 (]R 3 ), which 
implies the existence of a continuous embedding 

Wj [£«(R 3 )p, se[2,6] 

(see Edmunds and Evans [B, Theorem V.3.7). We denote by H„ (O) the closure 
of [Co°(fi)] TO in The following propositions follow in a similar way to Lemma 
2.1 and Lemma 2.2 in || and help to satisfy our needs in Theorem |l|. 

Proposition 1. Let a v € L 2 loc (Q), v = 1,2,3, and W, |B| <E L 3 / 2 (ft). J7ien 
Ply = Po + W is defined as a form sum with form domain Ti\ (f2). 

Proposition 2. Let a, V(= W), |B| satisfy the conditions of Proposition [| and 
suppose that the a v are locally Lipschitz on fi. Then Dy is essentially self-adjoint 
on [C(f(fi)] 4 . If furthermore V £ L 3 (fi), then B v has domain H^ (fi). 

At the end of §2 of |J conditions, consistent with those in Proposition ||, are 
given which ensure that IBV has essential spectrum K \ (—1, 1). 

The assumptions made in Propositions [I] and ^ are in fact sufficient for there 
to exist only a finite number of negative eigenvalues of ¥w(ty (and in (—1, 1) for 
Dy(O)). To consider cases in which |B| is not in L 3 / 2 (f2) we make use of an estimate 
given by Shen in |2(]] for the trace Mi(P\y, R 3 )- The assumptions that this requires 
are sufficient to ensure that Pw is defined as a form sum. To be specific, let 

*„(x) :=sup|;>0:; 2 (ja J q{ () l^(y)l P ^ <1 
where Q(x, I) denotes the cube in R 3 center x and side I. Define 

This is Shen's "effective" magnetic field. We recall that the need to replace |B| 
by some screened version in Lieb-Thirring inequalities for Pauli operators with 
non-homogeneous magnetic fields was demonstrated by Erdos in H, and this had 
motivated Sobolev in ]2^ | to obtain Lieb-Thirring estimates similar in form to that 
obtained for constant fields in Lieb, Solovej, and Yngvason |l5) but with |B| re- 
placed by an "effective" magnetic field. Another result of this kind is established in 
j2| . The novelty of Shen's approach lies in the simple and natural way in which his 
bp is constructed. In [po| , Remark 1.4, the following result is proved for ft = R 3 , 
but the proof is valid in general. It also holds if b p is Sobolev's effective magnetic 
field, or that of Bugliaro et al. 

Proposition 3. Let a € L 2 oc (f2,R 3 ) and suppose that for any p > 3/2 and 7 > 1, 

W^ +3/2 , bl /2 Wl e L^Cl). Then Pvy is defined as a form sum with form domain 
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In the case when B has constant direction, there is another result of Shen [ p0| 
which fits our purpose, namely his estimate for the counting function 7V(Pyi/+A, ]R 3 ). 
Let B(x) = (0, 0, B(x)), where x = {x,xz),x € M 2 , and define 

l p (x):=sup^l>0:l 2 ^\B(v)\>etoj <1 

where S(x, I) denotes the square in R 2 center x and side /. Define 

1 



(10) W ' MxW 

The following result is proved in [^0|, Remark 1.4]. 

Proposition 4. Let a € Lf oc (£l,M. 3 ) and suppose that for any p > 1 and 7 > 1/2, 
Wl +3/2 , b p WZ +1/2 £ L^O). Then, with B = (0,0,5), F w is defined as a form 
sum with form domain 7i\ . 

In the rest of this section we present preparatory results for subsequent sections. 
Our initial assumptions are as follows: 

(^4i) Let Q = U^LjQfc, where {Qk Q ^ : k = 1, ...,K} is a finite collection of 
non-overlapping congruent cubes whose edges are parallel to the coordinate 
axes and of length r. The set Q is fixed, but the side lengths r and number 
K will depend on H in due course. 

(A2) Let Vo and Wq be piecewise constant functions, taking constant values in each 
Qk, and zero outside Q. 

(A 3 ) Assume that B is continuous. Define B° on Q by B°(x) = B(xfc), x e Qk, 
where x^ is the center of the cube Qk, and let B°(x) = for x ^ Q. Choosef] 
gauges a and a for B and B°, respectively, such that for every Qk 

(11) max |a(x) — a(x)| < Cra r , oy :— max |B(x) — B(y)|. 

xGQfc |x-y|<r 

Define 

(12) fB 7 (M,«) :=/3 7 f 6(x) 

Jn 

for 7 > where 

1 - 1 



w_(x)^H2^[2tt(x)+w(x)]l +l 



k=l 



f/x 



Pi ■= ^2 y o (1 - 1)-'* 7 *. 

For simplicity, we shall write 23(6, v, £1) for S$o(&, v, O) when there is no danger of 
confusion. We are interested in the asymptotic behavior of the Riesz means 

k 

where each Afc is a negative eigenvalue of the Pauli operator Pvk(£^), or each A& is 
an eigenvalue of the Dirac operator in (—1, 1). 

We make frequent use of Proposition 3.2 of Sobolev which follows from 
results of Colin de Verdiere ||. It states that there exists a constant C such that 
for any 5 e (0, |) 

(13) h 3 N(P (B°) + A, Q k ) < 58(M|B°|, A, Q k ) 



1 Scc the Appendix for details. 
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and 

Cft 2 

(14) h 3 N(F (B°) + A, Q k ) > (1 - S) 3( B(fih\B°\, A + — , Q fc ). 

or 

Here Po(B°)(Qfc) is the Dirichlet operator on Qk, and iV(Po(Bo) + A, Qk) the num- 
ber of its negative eigenvalues below —A < 0. 

Lemma 1. Define a°° := a a. For a// € [C$°(Qk)] 4 and 6 > 



,, A Jv (B) - l/ ) 2 - 1 + Wo - (Po(B°) + W )h]<f>, 4>) Qk 
[L0 > < e(P o (B°)I 2 0,0) Qfc + M 2 (1 + |)(|a°°| 2 0,0) Qfc . 

Proo/. Set D a := fV — and D Q = D&. We have that 

(py„-y ) 2 -i + ^, 



= ([(a- J D a +/3) 2 -l + W ] 
= ||[a-( J Do-^ oo )^ll 2 3fc + (^o0»Q fc 

= ||a • D o m k - 2vMt[([a ■ D ]cj>, [a ■ a°°]0) Qfc ] + M 2 || [a • «°°]0||^ fc + (W </>, 
= (P Wo (B°)I 2 0, 0) Qfc - 2/ifRe[([a • [a ■ a°°]<£) Qfc ] + M 2 || [a • a°°W Qk 

since ||ac - Do^HL = (P (B°)I 2 </>, 0) Qfc . Therefore, 



| ([(D Vo - V ) 2 - 1 + Wo - P Wo (B°)I 2 ]0, 0) ( 
< 6(F (B°)h<t>,<j>) + (1 + £)// 2 ||[a • a 001 ^ 112 



» 

which completes the proof. □ 



It follows from Lemma @ and (0) that for e [C5°(Q fe )] 4 and 9 e (0, 1) 

([(©y o (B)-y o ) 2 -l + ^ O ]0»Qfc < 

1 Dj (1 + (S>)(P (B°)I a 6 ^) Qfc + W \\<f>\\ 2 Qk + ^»W r U\\l k 

and 

([(©Vo(B)-V o ) 2 -l + VK o ]0,0) Qfc > 
{ '> (1 - 0)(P o (B<%0, ^) Qfc + WoU\\% k - ^W r Uf Qk . 

Hence, from (JX3j) and (113) we have that 



(18) 
and 
(19) 



ift 3 iV([(B yo (B) - Vo) 2 - 1 + Wo]; Qfc) > 

(1 - 8) 3 <8(»h\B°\, ^ + j$^fr 2 al + $$;Q k ) 



±h 6 N([(B Vo (B)-V ) 2 -l + W ];Q k ) < 
^h\B%^ - j^^alQk). 

Similar estimates follow for the Pauli operator: 

r9n , h 3 N(F Wo (B);Q k )> 

( ' (1 - 5?<B(»h\B% ^ + ^M 2 r 2 a, 2 + Q fc ) 

and 

, on h 3 N(P Wo (B);Q k ) < 



PAULI AND DIRAC OPERATORS WITH LARGE MAGNETIC FIELDS 



7 



To estimate error terms, we use the following inequalities for 03 = Q3o established 
in Sobolev : for any A > and any subset G of £1 



(22) 

and 
(23) 
where 



\<B(b,U 1 +\ > G)-<B(b,U 2 + \,G)\ 
<M(b, Ui - U 2 ,G) + J G \Ux - C/ 2 | 1/2 (|C/i| + |C/ 2 |)rfx 

<M(b, U!-U 2 , G) +N{U 1 - U 2 , G)h{Jf(Ui -U 2 ,G)i + JV(U 2 , G)l}, 



\<B(bi,U,G)-<B(b 2 ,U,G)\< 

M(\h -b 2 \,U,G) + M(\h - b 2 \, U, G)ijV(U, G)i 
+M(\bi - b 2 \,U, G)iM{U, G)i + M(\h - b 2 \,U, G)iM{U, G)* 



M{b,U,G):= / b(x)\U(x)\*dx, b > 0, Af(U,G) := / \U(x)\*dx. 

JG JG 



Lemma 2. If r — Ah, then 
(24) 

IQfcl- 1 5B(M|B°|, ^ + G!^g^ + G 2 ^, Q fc ) - B(/ift|B|, W , Qk) 



i+0 

A 2 ,,2fc2^.2 



<^h\B°\{8 + ^^h?a 2 r -r W/ 
Proof. With U 2 :=W and 



—p?h?al + ^4532)2 + (^ftoy)i + fj,ha r . 



Ui := + C.^l + C 2 K 



and 

Therefore, 



1 + 0) ^r 2 ' 

f 

^-H 2 h 2 al + -jfat -. F(9, A, fih, S, a r ) 

x e Qk- 



1 



< 



IB - B°| < a r 



M{[ih\¥, - B°|, W , Q k )<nha r \Q k 



M^h\B°\,U 1 - U 2 , Q fe )< M ft|B°|F(0, A, /i, ft, <5, oy) 2 |Q fc |, M(U 2 ,Q k )<\Q k \, 
and 

A^df/i - t^l, Qfc)<P(0, -4, M> S,a r )i \Q k \. 
The lemma follows from (|2C|) and (EHI). □ 



3. The Pauli operator: \ih < constant. 

We assume that (Ai)-(A 3 ) hold in §2 with V = 0. It follows that F Wo (fl) < 
)^ 1 PvK (Qfc) m the form sense (see [[|),§XI.2.2]), which implies that 



A" 



(25) 



N(F Wo ,Sl)>y N(F 



k=\ 
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and, since 

M 1 (¥ Wo + A, O) = - f Q °° f>dN(P Wo +X + t,Q) 

(26) 

= 7 / °° V- x N(W Wa + A + t, Q)dt, 
for all A > 0, we have for all 7 > 

K 

(27) M y (¥ Wa +\ 7 il)>Y / M 1 {¥ Wo + A, Q k ). 

k=l 

Note that (§|), and hence @, hold for A = if 2V(Pw„, Q) < 00, which will be the 
case in Theorem [j] below, but not in Theorem ^ and Theorem [| In Theorem || and 
Theorem ^ all we know is that N(Pw + A, Q) = 0(A _1 ) and 0(A~ 2 ) respectively, 
and as a consequence, we can claim that ( p6| ) holds for A = only if 7 > 1 in 
the first case and 7 > ^ in the second. We are mainly interested in the remaining 
values of 7 in Theorems || and ^, and take A > 0. 

Also, note that P\y > — (Wo)- > — A, say, and the inequalities ( ^p| ) and (21) 
with Wq replaced by Wq + A hold uniformly for A E [0, A]. Moreover, 



/>oc 

(28) B 7 (&,W,G) =7 / f r - 1 Bo(6,W + t,G)dt 

Hence, from (18), (|2^), and (26), we see that when p,h<l 
h 3 M^(F Wo + A;fi) > (l-£) 3 *8 7 (^|B|,T^o + A,Q) 



-0 



{((? + ^ + W) 4 + (* + TT^ + A^) 1 + + Tr} |Q|) • 



On allowing 7i — > 0, — > 0, A — > 00, and S — > in that order we have that for all 
7 > 

(29) liminf{fi 3 M 7 (Pvy + A; O) - QS 7 (/zfi|B|, Wq + A, fi)} > 0. 

h— *o 

To prove the reverse inequality in ( |29| ) , we proceed in the manner of Y. Colin de 
Verdiere || and Sobolev Let the interior of each Qk be denoted by int(Qk)- 
Set 

(30) S :=M. 3 \U%_ 1 int(Q k ), S pr := {1 £ I 3 : dist(x. S) < pr} 

for some p £ (0,1). Construct a partition of unity {V'fclfcLo subordinate to the 
covering U^ =1 int(Qk) U S pr of R 3 : 

(i) ^0 € C°°(R 3 ), ^eC °°(mt(Q fc )),A;>l, 

(31) (**) ELoV'N 1 . 

Ef=o l v V'A(af)| 2 <(/w) _2 Xpr ) 

where % pr is the characteristic function of Q pr :— {x € M 3 : dist(x, Q) < pr}. Note 
that 

(32) \S pr nQ\< P \Q\, \Q pr \Q\<pr. 
Then, for every / 6 [Cq 30 ^)] 2 , a calculation yields 

K 

(Polto/, ifo/) + ]T( P ^< ^ ( P °/' /) + Cft 2 (pr)- 2 ( Xpr /, /) 



fe=X 
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and, as a consequence, 
(33) 



(Pwo(n)/,/) > ([Pw (s P r no) - ch 2 (pr)- 2 X pr}i>of,M) 



By the minimax principle, A„(V>fcPw (Q)V'fc) > A n (Pwn(Q)) for fc = 0,1,..., if, 
(see (21) of pi). Therefore, it follows from (pi and (||) that for all 7 > 



M 7 (P Wo + A, ft) < M 7 (P Wo + A - Cft 2 (pr)- 2 x P ,., S pr H 0) 
1 j +Ef=i^7(^o+A-^ 2 (pr)- 2 Xpr ,Q fc ). 

From (plj), we have in the case that /j,h<l, and uniformly for t G [0, A], 

7i 3 AT(P Wo + A + t - Ch 2 {pr)- 2 ,Q k ) 

S -DU^"! 13 I) 1=0 9(1-8) V r °V (iJ9) p i r 'i 1 Wk) 

< B(^fi|B|, Wo + A + i, Q fc ) + o({(0 + ^a 2 + 3^)* 

+ (<? + + + (trr)i + ^}|Q fc |) 

as in Lemma ||. Hence, on using (|26|), 

^ Ef =1 M 7 (P Wo + A - C?i 2 (pr)- 2 , Q fc ) 

(35) < 58 7 (^ft|B|, Wo + A, Q) + o((fl + ^a 2 + 
+ (0+^ 2 + ^)l+a|+ ( 

For the first term on the right side of the inequality ( |34| ) , we use the Cwikel-Lieb- 
Rozenblum inequality for a Schrodinger operator with a magnetic fieloj^] to derive 
the estimate 

ft 3 A/ 7 (P^ + A - Ch 2 {pr)- 2 x P r, s pr n ft) 

< Js^nnWI -W + ^ Xpr ]J + ^x 

(36) < J^nodl^ol + |B|]I+^ + [A] 2 7+3 Xpr)d: 



x 



where 



<fn\ Q |B|i+^x + / SprnQ [|W | + |B|]l+^x+ [^y] 2 ^ 3 |Q pr n S pr \ 
<fn\Q IBI^dx + K{pr) + (Ap) l 7+3 [p+pr] 

K{pr)= [ [|W o | + |B|]i+^x->0, as p 0, 



's pr .nc 

uniformly in fi, and A, in view of (|32|). From (|34|), (|35|), and (|36"1), on allowing h — > 0, 
$ — ► 0, yl — ► 00, and p — * 0, in that order, it follows that 

(37) limsup{n 3 M 7 (P Wo + A, ft) - Q5 7 (jufi|B|, W + A,ft)} < / |B(x)|2+T(fa;. 

Since Wo = outside Q, then Q can be chosen such that the right-hand side of 
@ is arbitrarily small if |B| e L 3 / 2+7 (ft). From this fact and @, we conclude 
that for all 7 > 

(38) lim {ft 3 M 7 (Pvy + A; ft) - <B 7 (/i/i|B|, Wo + A, ft)} = 

Now, we are in a position to prove the following general result for a Pauli operator 
with an electric potential that is not assumed to be piecewise constant. 



2 Sec Theorem 2.15 of Avron et al. nil or Lemma 4.1 of 
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Theorem 1. Suppose that 

(i) B is continuous on Q, 

(ii) w, |B| e Li+^n), 

(iii) [ih < constant. 
Then, for all 7 > 0, 

lxm{H 3 MJW> w ,Q) - 35 7 (/ifi|B|, = 

h— >0 

Proof. Given e > 0, there exists a finite collection of non-overlapping congruent 
cubes {Q k }k=i for which (^lM^a) of §2 hold (for V = 0) with 

(39) \\W~W Q \\s + ^ n <e, ||B-B ||3 +7 , n <e. 

Note that Wo depends on e and hence so does the lower bound A of Pw - Let 
rj G (0, 1) and set 

P w = T x +T 2 , 

(40) T x := (1 - r/)P + W - rififux ■ B° 
T 2 := r^Po + Woe + wha ■ B° 

for Woo ■= W - W Q . Then, 

Tx= [a ■ (-ihV - jla)] 2 + W - rj^itkr ■ B° 
T 2 = v H (B) + W QC -wha-B oc 

where 

h:=(l-t))*h, ft := (1 -T])*ii, and B°° := B - B°. 

We have that 

Tx > Sx := [a ■ {-ihV - fm)] 2 + W - n^h\B°\ 

and 

T 2 > S 2 := nH Q {B) + W^ - W h\B°°\. 

It then follows that 

(41) M 7 (P w , n) < M 7 (S 1 ! , n) + M 7 (5a , fi) . 

From ( |38| ) with A = (which, according to the remark after (|27|), is allowed since 

N(Sx,n) < 00) 

(42) lim{/i 3 M 7 (Si,n) -<8 7 (£/j|B|,Wo - Jtfifi|B°|, fi)} = 0. 
and by Theorem 2.15 of Avron et al. Q 

(43) fi 3 M 7 (S 2 ,fi)< / [77- 'Woo -^|B 0O |] i+7 dx<( e /7 ? )l+^. 

~ Jo 

Let 7 = 0. On using (|22]), we have that 

|tB(M|B|, Wo ~ 77^ft|B |, O) - «8(/ift|B| s W,n)| 

(44) ^/nlBlNWool+^B^ds 

+ ||W 0C +r ? |B ||||^{||T^ 0O + r ? |B°||| | ^ + ||W||3,J. 

It follows from @-@ that 

limsup{n 3 iV(Pvi/;^)- < B(^|B|,T^,f7)} < 0. 
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The cases 7 > follow similarly from ( p6| ) and the inequalities (2.23) - (2.25) for B 7 
in [ p3[ . The reverse inequality is obtained by choosing r\ S (—1,0) and repeating 
the argument. □ 

If |B| is not in L 3 / 2 (fi), there may be an infinite number of negative eigenvalues. 
In this case we have the following 

Theorem 2. Suppose that 

1. B is continuous, 

2. |B|, w e L°°{n), 

3. for p > 3/2, W 5 / 2 , b 3 p /2 W e L 1 ^), where b p is defined in (gj, 

4. ph < constant. 

Then, for all 7 € [0, 1) and A > 0, 

(45) lim{ft 3 Af 7 (P M / + A;fi)- ( 8 7 (/ift|B|,M/ + A,fi)} =0. 
If 7 > 1 j (fPD wii/i A = is proved in p3| . 

Proof. We first note that Pjy is properly defined as a form sum by Proposition |[ 
Also, from Theorem 1.1] for A > 0, 

7V(P W + A,K 3 ) < lAr^^R 3 ) 

(46) < IX]- 1 ^- 3 [ w! /2 dx + ,?/ 2 /r 3 / 2 [ 6 3 / 2 W_dx} 
and hence, for any 7 > and A > 

(47) M 7 (P W + A,f7)<|A|- 1 {fi- 3 / iy! /2 dx + M 3 / 2 ft- 3 / 2 / 6 3/2 iy_dx} 

Jo Jo 

Suppose that W , B°, Q = U^ x Q fc satisfy (Ai)-(A 3 ) of §2. Then (j29j) follows as 
before, and so do ( |34| ) and (|35| ), with Wo replaced by Wq + A. The remaining term 
on the right-hand side of (|34j) is estimated by (fl6|). Since (Wo + A)_ < (Wo)-, we 
have, with r = Ah, 

\x\h 3 M 7 (p Wa + a - ch 2 (pr)- 2 Xp r, s pr n n) 

< Js P ,.no ([^0 - Cft 2 (pr)- 2 Xpr ]i /2 + b 3 p /2 [W - Ch 2 (pr)- 2 Xf >r}-)dx 

</s P ,.no([^o] 5 _ /2 + 6r[^o]-)rfx 

+ (ft/pr) 5 |S pr n Q pr \ + (h/pr) 2 J SprnQpr b 3 p /2 dx 
<K(pr) + {l/Apf{p + pr) + (l/Ap) 2 (p + pr), 



where K(pr) — > as p — > 0, uniformly in h and A, by (|32|). It follows as for (38) 
that 

(48) lim {h z M 1 {V Wo + A; fi) - <8 7 (/xft|B|, W + A, Q)} = 0. 

For a general W we proceed in a similar way to Sobolev in p3[ |, using Shen's 
estimate (|4^). 

Let e > and choose non-overlapping cubes Q^, k = X, K, with sides parallel 
to the co-ordinate planes, and a piecewise constant Wo which is constant in each 
Qk, is zero outside Q = Uj^Qfc and 

(49) / \W -W \ 3/2 dx<e, I b p \W - W | 1/2 dx < e. 
Jn Jn 
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Note that assumptions 2. and 3. imply that W 3/2 , b p W 1/2 £ L x (f2) and 

(50) / \W -W \ 5/2 dx<e, [ b 3/2 \W - W \dx<e. 
Jn Jn 

We prove the result for 7 = 0, the proof for 7 > being similar. Set N\(W) 
N(¥w + A, ^)- Then, we have the Weyl inequality 

N X (W) > N X (W /(1 + 0) - N X {[W ~ Wo]/C) 
for any ( E (0, 1). By © and © 

(51) h 3 N x ([W-W }/()< W-^eC-W + eC 1 ). 
Also, from @, 

lim {h 3 N x (W /(l + 0) - «8 (M|B|, W /(l + C) + A, fi)} = 



and from (|22j) we have 

|B (Mft|B|, W /(l + + A, 0) - fB (M|B|, W + A, f!)|<C 1/2 + e 1/2 . 

Hence, 

liminf^o {h 3 N x (W) - Q3 (M|B|, W + A, SI) } > 
-lAl- 1 ^- 5 / 2 + er 1 ) - C 1/2 - £ 1/2 - 
Let e — > 0, C — > in that order to get 

liminf {h 3 N x (W) - ^B (fj,h\B\,W + A, fi)} > 0. 
The proof of the upper bound is similar. We use 

N X (W)< N x (W /(l-())+N x ([W- W ]/C) 
and argue as before. □ 

The following result can be established similarly. 

Theorem 3. Suppose that 

(i) W , B are continuous, 

(ii) |W(x)| — * uniformly as |x| — > 00 in O, and 

(iii) fih < constant. 

Then, ZioZds /or a// 7 £ [0, 1) and A > 0. 

Proof. The operator Pvk(^) is properly defined as a self-adjoint operator on the 
domain of Po(fi) since it is a bounded perturbation of Po(f2). The point to note 
is that [W + A]_ is compactly supported in ft. On choosing Q in (Ax) of §2 
to contain this support, and the piecewise constant function Wo to be such that 
Woo = W - Wo > in Q, we have for 



w 



+ 2A = {(1 - ry)P + W + A} + {nP + + A} 



that R 2 > for all 77 € (0, 1). The proof follows easily from fl48|). □ 
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4. The DlRAC operator: ah < constant 

Assume (Ai)-(A 3 ) of §2 with W = 0. Set N(S,ft,I) := #{A„(S) : X n (S) € I}, 
the number of eigenvalues of the operator S(f2) in the interval I, and M T (5; fi, I) — 
Ea„(S)£/ l^™( 5 ')l 7 - Then, with V > in Q fc , 

7V(© Vo , Q fc , (-1, 1)) + JV(D_ Vo , Q fe) (-1, 1)) 

= N(B V(n Q k , (-1, 1)) + iV(Dy , Q fc , (-1 + 2V , 1 + 2V )) 

= N(Bv ,Q k ,(-l,l + 2V ))+E, 



where 



iV(Dy o ,Q fc ,(-l + 2F ,l)) = if Vb < 1, 
-JV(Dv ,Qfc,[l,-l + 2Vo)) if V& > 1. 



For M := - 1 = P0I2 and A° := (1 + V ) 2 - 1, 

AT(IV , Q fc , (-1, 1 + 2V )) = N{Rq, Q k , [-1, A°)), 

which implies that for any 77 > 

\N(D Vo ,Q k , (-1, 1)) + iV(B_ FoJ Q fc , (-1, 1)) - 7V(R , Q fc) [-1, A?))| 
1 ' <N(Ro,Qk,[-M-i + r))) 

where A^ := (1 — Vb) 2 — 1. Inequality ( |52| ) holds as well for Vb < with 

(53) A° := (1 + |Vb|) 2 - 1, A° a := (1 - |V |) 2 - 1. 

For R Vo = Ry (B) := Dy 2 - 1, we have My (fi) < ®k =1 M Vo (Q k ). Since 

(54) iV(D Vo ,O,Z)=iV(My ,O)=iV(My ,O,(-oo,0)), I =(-1,1), 

we may apply the minimax principle to Ry in order to estimate iV(By- D , fi, I). In 
particular, 

7V( D y , 0, 7) + AT(D_ Fo , n, J) > Ek=i{N(Rv a ,Qk) + N(R- Va , Q fc )} 

= J2k=i{N(Dv ,Q k ,I) + N(B- Vo ,Q k ,I)}. 

Hence, from @, @, and (H) 

in 3 {iv(Byo , n, /) + ^(©- v , O, I)} 

> l ft3 Ef=i{^(Ko, 0*, [-1, A?)) - iV(JRo, Qk, [-1, XU + V))} 

> (1 - <5) 3 <B Wl, ^ + + Q) 
-5B(^|B°|,^^-^S,Q) 

= (1 - 5) 3 5B(M|B|, -A?, Q) - 95(M|B|, -AL X , Q) 
_0({(e + rj+^ + ^)i + (<? + »?+ ^ + ^)i+^+a r }|Q|) 

for r = Aft, on using Lemma [T^ with nh<l. It follows as for (29) that 



(55) 



lim inf ft ^o {5ft 3 [iV(Dvb ,0,7)+ iV(B_ M , , n, I)] 



-«8(Aifi|B|, -A?, fi) + Q5(/xft|B|, -A° !, O)} > 0. 
Note that X°_ 1 < if \Vq\ < 2, and consequently, 
(56) 5B(/x?i|B|,-A° 1 ,Q fc )=Q if |V | < 2. 



14 A. A. BALINSKY, W. D. EVANS, AND R. T. LEWIS 

To establish the reverse inequality in ([55]), we first note that for {4>k\k=o gi ven hi 
© and all /€ [C °°(f!)] 4 

h 3 

Since EfeLo ^fe = 1> ^ follows that 

= Ovo 2 / + ^ 2 Ef=olV^| 2 /- 

Therefore, by © 



(R % («)/,/) > ELi([Kv (Qfc) - j&WkfM) 

+([Rv (Spr nil)- £$x P r]M,M). 



(57) 

Let a 2 := Ch 2 /(pr) 2 , where C is the constant in (|57|). Then, 



iV(Ry - a 2 , Q fc ) - iV(© Vo , Q fc , (-Vl + a 2 , Vl + a 2 )), 
and on repeating the argument leading to (52), we have for any rj > 0, with A^ :1 (a) = 
(x/T+^ilVbl) 2 -!, 
(58) 

±h 3 {N(R Vo - a 2 , Q fc ) + N(R_ Va - a 2 , Q fe )} 

< ±h 3 {N(R Q - A?(a), Q fe ) + 7V(K - A° t (a) - tj, Q fc )} 

< X(»H\B°\, -^M _ , Qfe) + <B(^|B°|, _ _ (CW , Qfc) 

< 95(/ifi|B|,-Af, Qfe) + «8(Mft|B|, -A° 1; Q fe ) 

+ (# + ?? + ^ + ^+) § + +^r}|Qfe|) 

by ( |Tg| ) and (p4|), since r = AS. From Lemma 4.1 of Q 
(59) 

^ 3 iV(Ky - ffixmSer n fi )< Js pr nJ y °l S + l^o| 3 + |B|I + j^x P r]d^ 

< In\Q |B(x)|ldx + K(pr) + {Ap)- 3 [p + pr] 

where 

K(pr):= f [|Vb|5 + |Vb| 3 |B|i]dx^0, as p -> 0, 

Js pr nQ 

uniformly in ft and A by @. As in on using Q, it follows from @, (||) 
and @ that 

limsup^ |ift 3 [iV(© yQ ,fi,/) + iV(D_ yo ,r!,/)] 
k ^ -®(Atfi|B|, -A?, fl) - ®(Atfi|B|, -A° !, fi)} < 

for J = (-1,1). 

Now, we are in a position to prove the first theorem in this section. 

Theorem 4. Suppose that 

(i) B is continuous on Q, 

(ii) V,^ 2 ,|B|eLl+^) ; 

(iii) |{x G fi : |7(x)| > 2}| = 0, 

(iv) ph < constant. 
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Then, with I v := (—1 + rj, 1 — rf), rj G (0, 1), and any 7 £ [0, 1], 

lim \im{lh 3 [M. y {B v -n,I v ) + M 7 (B>_y; fi, /„)] - 58 7 (>?j|B|, -Ai, O)} = 

77— >0 ?i^0 2 

w/iere Ai := 2|V| + V 2 . 

Proof. We prove the result for 7 = 0, the other cases being similar. As in the proof 
of Theorem |l|, given e > choose a piecewise constant function Vo taking constant 
values in non-overlapping congruent cubes Qk, k = 1, . . . , K, and such that 

(61) \\V-V \\ hQ <e, \\V - V \\ 3 ,a < e ||B - B°|| f ^ < e 

in which Vo = 0, B° = outside Q := U% =1 Q k . Thus, §|) and © hold for Vo, 
and, indeed, for any function taking constant values in each Qk and zero outside 
Q. Let 77 e (0, 1) and set 

By = (1 - n){a ■ (f V - /xa)} + P + V + v{a ■ (f V - /ia)} 

(62) = {a.(|v-Aa)} + /3 + Vb} + {r 7 a-(fV-/ia) + y oo } 
= : T!+T 2 

where ft := (1 - rj)K, fj, := (1 - r?)^, and Voo = V - V . Then, for e [C^°(fi)] 4 , 

(63) ([B y - 1]0,0) = IIT^H 2 + 2SU(T 1 <^ ) T 2 <f>) + \\T 2 ^\\ 2 - U\\ 2 . 

With D a := jV - jla and D a := |V — ya, we have 

(T!0, T 2 0) = ({a ■ £> a + 13 + V }<P, {^-a ■ D a + V^) 
= ^-WT^ + iT^iV^-j^iP + VoM. 

Since 9tc(a ■ £>„(/>, 04>) = 0, 

2£Re(Ti&T 2 = 

^WT^f + 23te[(T 1 & t^) - ^(T^, Vb^) - ^((/3 + 

Hence, 

|2<Ke(T^,T 2 0) - t^IIT^II 2 + ^||0|| 2 | 
'" 1 < 20||T 1 0|| 2 + HWV^f + ^HVo^ll 2 ) + t^(| 

Also, 

||T 2 0|| 2 = 77 2 ||a • Da^W 2 + 2 V mt{a ■ D*<j>, V^cp) + || V^f 
which implies that for any 6 > 

^ < ffy 2 ([#0(B) - jtzTlCT • B]I 2 0, 0) + (1 + i)|^oo0|| 2 . 

On substituting (^) and fl65| ) in ([33]), we have that 

([D 2 , - 1]0, 4>)> (1 + i 2 ^ - 20)||T 1 0|| 2 - i(ll^0H 2 + (T^rll W) 
-j^(|F |^ + 77 2 (1 - 0)([ff o (B) - /ifta • B]I a ^, 0) 

-(l + i)||Ko0|| 2 -||0|| 2 - T^ll^ll 2 

( 66 ) > (l + ^-28)(lB 2 Vo -l}t^)- I ^\\V 0( j ) \\ 2 

-i*-(\V \^ 4>) - r? 2 (l - 0) M ft(|B°|<M) 
+r? 2 (l - 6)([H (B) - ^|B°°|]1 2 0, 0) 
-(l + f)(^>,0)-20||0|| 2 
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where Dy = T i = a-(|V-/ia) + /3+Fo. Whence, on choosing 9 = 77/(1-77) < 1/2, 
mfi_ _ f-i _ jn. 



"V 



(1 - 1^) > K 1 - W) 2 + 2|Vb|) - V 2 ^\B°\ 



( 67 ) +\[H (B)-txh\B°°\]h-^V2 



with 



Si = Si(V ) := D 2 / - 1 - $, $ := -^-(^ + 2|V |) + 77 2 ^|B°| 
" 1 — 



and 

S 2 := ^[ifo(B) - ^|B-| - ^^]I 2 . 
2 77^ 

Since is piecewise constant and vanishes outside Q, the analysis leading to in- 
equality (58) holds with Ry = Dy g — 1 and a 2 = $, and this yields the estimate 

p 3 {jV(5i(F ), O) + iV(Si(-7 ), n)} 
< 5B(M|B|, -A?, fi) + ®(M|B|, -A° l5 n) 

+O([0 + 77^ + + [0 + V* + + Or + ff P ) 

where A^ = (1 ± |Vb|) 2 - 1, 9 is arbitrary and r = Ah. Also, by (§|) 
|93(M|B|,-A° ±1 ,fi)-£(M|B|,-A ±1 ,f!)| = 0(11^111 + ||V M || 3i0 + ||B|U n )<e 1 / 2 

with A±i := (1 ± |F|) 2 - 1. It follows that 

Umsup^ {|h 3 [iV(5i(yo) ) 0) + iV(5 1 (-Vo),0)] 

-Q5(M|B|,-A 1 ,Q)-5B(A/i|B|,-A_ 1 ,fi)}<e5+r 7 i. 

A similar argument yields the reverse inequality (c.f.(|55|)) 

liminf r ^o{^ 3 {A r (^i(K ) ),^) + A r (5i(-Vb),r!)} 

(69) -B(M|B|,-Ai,fi) + B(M|B|,-A_i,fi)}>- (e* + 7? j). 

Note that A_i = (1 - |V|) 2 - 1 < for |V| < 2, and also [2kfih\B\ - A_ x ]_ < V 2 , 
and is zero for 2fe/i/i|B| > A_i. Therefore, we have 

B(£ft|B|,-A_i,n)< / |V(a:)|[|B(a;)| + |V(a;)|]da; = 

from assumption (iii). Hence, as H —* 0, 

(70) ^ 3 {A^(5i(Yo), 0) + iV(5i(-F ), O)} - B(M|B|, -A X) fi) 
= 0(e^+77^). 

Finally, by ||, Lemma 4.1, 

(71) /7 3 A^(S* 2 ,fi)< / [|B°°|i +?7-5V^]dx<ei +e 3 ?7-i 
It follows from @ ,@, and @ that, with J,, = (-1 + 77, 1 - 77), 

hmsup{-/2 3 [iV(Bv,fi ) J^) + A/'(irj)_ v ,,fi,^)] - «8(/ifi|B|,-Ai,fi)}<?7*. 

The reverse inequality is proved similarly, with 77 < in (|62]). The proof is complete. 

□ 
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The next theorem is the analogue of Theorem [2] in §3. 

Theorem 5. Suppose that 

1. B is continuous, 

2. |B|, V€L°°(n), 

3. for some p > 3/2, V 5 , b^V 2 € L^O), 

4. \xh < constant, 

5. |{x e n : |F(x)| > 2} I = 0. 
Then, for all 7 G [0, 1], 

h 3 

lim lim {— \MJB v ,n,I„) + MJB- V ,tl,I v )] - 5B~(^|B|, -Ai,ft)} =0, 

77— >0 R— >0 2 

where I v = (-1 + 77, 1 - 77), 77 € (0,1) and Ai =2\V\ + V 2 . 

Proof. The only difference with the proof of Theorem ^ is in the way the error 
terms involving S pr n are estimated to give the analogue of . These are now 
dependent on (^) rather than on the magnetic CLR inequality. We again prove it 
for 7 = 0. From (By - V) 2 = (P + l)I 2j it readily follows that, for any C > 0, 

and, for any A < 1, on choosing £ = A/(4 — A), 

B v 2 - 1 + A > (3/4){P + A - (4/A)^ 2 }. 

Consequently, 

»Vo + A - ^(^iXpr- > (3/4) |P + A - (4/A)V 2 - C-^ Xpr 
and, from p^), with r = Ah, 

h 3 N(R Vo + a - c^xpr, 5^ n fi) 

< Is pr nn H V o + (l/Ap) 2 x P r} 5/2 + bl /2 [V 2 + l/(Ap) 2 Xpr }}dx 
<K(pr) + {l/Apf{p + pr) + [l/Ap) 2 (p + pr), 

as in the proof of Theorem |2| where K(pr) — » as p — > 0, uniformly in ft and A. 
This gives (^0|) as in Theorem 0. From (p6| ) we have, 

By 2 - 1 + 2A + 26 > (Si + A) + (S* 3 + A), 

where, with 8 = 77/(1 — 77), 

S 3 = V (l-2r,)F (B°°)-(l/e)V£. 

The Si + A term is handled as in the proof of Theorem || and the first part above. 
The S3 term on S pr n fl is again estimated by ( pj| ) . □ 

5. The Pauli Operator: magnetic fields with constant direction. 

First, we summarize some facts from |2^| that we will need below. Let B(x) — 
(0, 0, B(x\, X2)), x = (x, X3), x — (xi,X2)- Assume that B = B(x\, X2) is periodic, 
i.e., for some T\,T% 

B(x\,X2) — B(x\ + Ti, X2) — B(xi,X2 + T2). 
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For T 3 > denote by D (3) = xjL-JO, T k ) the fundamental domain of the lattice Y 
with vertices (Timi, T 2 m 2l T3?TJ3), mj E Z, j = 1, 2, 3. Choose the following vector 
potential a corresponding to B: for <f> = <j>{x\, X2) 

a(x) = {-d 2 <p,d 1 <j>,0), A<j> = B. 

Set 

B o= tTt^ Jo 1 Jo* B(x 1 ,x 2 )dx 1 dx 2 , B := (0, 0, B Q ), 
B\ = B — B , 

0o = ^-{axl+pxl), a + /? = 2, 

0= 00+01, 

where A0i = S x and 0i is periodic on D< 2 ) = [0,Ti) x [0,T 2 ). Note that the 
conditions above imply that B\ is periodic, / Q Tl J^ 2 B\dx\dx 2 = 0, and A0 O = -B . 
The /?ux of B across 

£)(2) 

IS 

$ = -!- / B(x 1 ,x 2 )dx 1 dx 2 . 

27T ,/£)(2) 

The physics of our system must be invariant under translations (xi,x 2 ) — > (xi + 
Ti, x 2 ), (xi, x 2 ) — > (xi,x 2 +T2), since the magnetic field is invariant under these 
translations. The magnetic potential a is not invariant, but since a and the trans- 
lated a, SLTj ,j = 1,2, say, give rise to the same magnetic field, there must be gauge 
transformations 

(r^)(x) :=e^ x V 3 .(x), j = 1,2, 

in which 

i*Ti(x) = u(xi +Ti,x 2 ), ur 2 (x) = u(xi,x 2 +T 2 ), 
which are unitary maps on the Hilbcrt space and 

(|v - im)(tju) = e if * (|v - fia Tj )u Tj , j = 1, 2. 

From the discussions above, we see that 

a = a°+a\ a fe := (-d 2 fe , 9i0 fe , 0), fc = 0, 1, 

where a 1 is periodic and 

a°(xi +Ti,x 2 ) = a°(x u x 2 ) + ^aTi(0, 1, 0) 
a°(xi,x 2 +T 2 ) = a°(xi, x 2 ) - ^/3T 2 (1, 0, 0). 

The choice /i(x) = -i ^§f^- x 2 and / 2 (x) = i ^° T2 Xl yield the "magnetic trans- 
lations" given by 

(t 1 u)(xi,x 2 ,x 3 ) = u(xi + T 1 ,x 2 ,x 3 )exp(-i^iaB Tix 2 /2h), 
(t 2 u)(x\,x 2 ,x 3 ) = u(xi, x 2 + T 2 , x 3 ) exp(i/j,pB T 2 xi/2h). 

They commute with the expressions Uk = —ihdk — [Mik, k = 1,2,3, and Q± := 
III ± »n 2 . If the flux of B satisfies the condition 

(72) ^/r 1 * = N G Z 

then Ti and r 2 commute and we can reduce our problem to one on the torus = 
M 3 /r. In this case, we define the operators IT^ = —ihdk — [tat, and Q± in L 2 (X^) 
with domains consisting of functions u € C°°(tt 3 ) which satisfy 

Tku^u, k = l,2, u(xi,x 2 ,X3 + T 3 ) = u(xi,x 2 ,x 3 ). 
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and denote them by Q±(X^), Ilk(X^). It is proved in pi| that these operators 
are closable, that each Hk is symmetric and Q± C Q^. Let 

A±(X^) := H (X^)) T ^hB 

= Q* ± (X^)Q ± {X^)+I1 2 (X^) 

where Hq = Hq(B), B = V x a, as defined in §1. We denote the closures of the 
operators by the same notation. The Pauli operator on the torus X^ 3 \ P(X' 3 )), is 
now defined by 

A+(X&) \ 
A^(X^) J 

and P W {X^) = Fq(X^) + W. We remind the reader that F W (Q),F W (Q) will 
always stand for the operator Pw with Dirichlet boundary conditions. 

Using a result^] due to Dubrovin and Novikov Q that if ±N > (N as defined 
in (|72])) A = is an eigenvalue of A±(X^) of multiplicity ±iV, Sobolev proves 

Lemma 3. (Sobolev p^j , Lemma 4.3, 7 = 0, d = 3) Let ffffij be satisfied and 
\B\ > k > 0. IfW is constant in and W- < 2[ihn, then 

\N(P W ,XW) - ^W1\N\\ < \N\. 

We shall follow Sobolev's strategy, which is again basically inspired by the 
method of Colin de Verdiere ||, but substantially modified to meet the needs of 
the constant direction magnetic field case, using the observations and results noted 
above. There will now be an infinity of negative eigenvalues in general, and so we 
can only expect results for N(Fw + A, O) with A > 0; in fact for M 1 iPw + A, fi), 
A > when 7 S [0,1/2]. When 7 > 1/2, M 7 (P w ,fi) is finite and its semi- 
classical limit is given by Sobolev in ^3|. Our analysis requires an estimate for 
h 3 N(F w + A, ft), A > 0, in which fi\B\ occurs linearly. A suitable result is obtained 
by Shen in |20| . In it N(Fw + A, fl) grows like A -1 / 2 as A — > 0+, and we are forced 
to take A > in our result for M 7 (Fw + A, f2) when 7 < 1/2 in view of (^6|). 
If fi is a bounded domain we prove in the next lemma that the A dependence in 
N(Pw + A, ft) is at most logarithmic; of course, this is still crude since the spec- 
trum of Pw(Q) is discrete in this case, but it does have interesting implications for 
Mj(Fw, ^) for 7 > 0, and, moreover, the technique used to prove it can be made 
to yield a similar estimate for unbounded domains; see Proposition]^ below. 

Lemma 4. Let Q = x [0,R), W- G L°°(Q), and B = (0,0,5) . Then, for 
A > and W := inf Q W{x) 

(73) N(P W + A, Q)<[(W- + l)e^- 1 log \~ 1 \]^,h^ 2 {W + A)i|Q| max \B(x)\. 

~ Q< 2 > 

Proof. We have that 

Pw = F^ + W- h 2 dl > P { q } +W- h 2 dl = P w 

(2) 

where Pq is the Pauli operator in two dimensions, namely 

3 See also Appendix A of Sobolev |23[ for a proof. 
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with 
and 



Af = Q*±Q± = H (2) TfihB 



fc=i 



the magnetic Schrodinger operator in R 2 . The eigenvalues of — fi, 2 <9f with Dirichlet 
boundary conditions on (0, R) are e m = (m7rS) 2 /i? 2 , m E N. It follows that the 



eigenvalues of are of the form A„(P n ; + W) + e m , where A„(P n ; + W) arc the 





eigenvalues of Pq 2 ^ + W" on (with Dirichlet boundary conditions). Hence, for 



(74) 



any AeK, 

N(P W + A,Q) < E™ ^( p o 2) + W + tm + A, Q^) 

= E m E - =1 ^(#o 2) + (-i) j 'm*b + W + e ro + a, q( 2 )). 

Since P*, 2) > 0, we have A^ = H [2) =F > 0. Hence, for any e E (0, 1) 

TnhB + W + e m + \> e(H {2) =p ^B) + M> + e m + A 
> e^ 2) - e^|B| + M> + e m + A. 

Replacing A± by something smaller seems rather wasteful, but will be seen to 
provide the correct scaling for fi\B\ by a suitable choice of e. It is reminiscent of a 
technique (running energy-scale rcnormalization) used in a paper of Lieb, Loss and 
Solovej Q. On substituting in (f74|), this yields 

(75) N{¥ w + A, Q) < 2 ^ A^(e^ 2) - €fih\B\ +W+e m + X, Q (2) ) 

m 

where the sum is over all m such that 

(76) e m = ^<(# + A)_. 

To estimate N{eH^ ] - efj,H\B\ -W-+ e m + A, Q^) we use the result of Rozenblum 
and Solomyak [[l9|, Theorem 2.4] with (in their notation) 

A = eH {2) + e m + A, B = -eh 2 A + e m + A, V = (efih\B\ + W-)xqm . 

Then, B E V as it generates a positivity preserving contractive semigroup with the 
(2, oo)-boundedness (ultracontractivity) property having kernel Q(t,x,y), where, 
on the diagonal, 

Q{t 7 x, x) — n g ^ ^ , 
2irn l! €t 

sec §2, for the terminology. The operator A is dominated by B in the sense 
that 

\e- tA ifj{x)\ < e' tts \ip(x)\ for a.e. x E K 2 , a: = (xi, a^), 

i.e., .A S 7T»(S) in the language of |l9), §2.3. It follows from [[§||, Theorem 2.4] 
with the choice G(z) = (z — k) + for some k > 0, that 

N(A -V)< jjL / °° f / 0W ^ e ~(-+A )t( ^ (x) _ fc)+dl 

(77) = 2OTV)^A7v(x) e " (£m+A) *(^)- fc )w 
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for <j)(x) := k/V(x) and 



(7M ,(1)= / (z-fe) + e-*^= /V^. 

o z h $ 



Now choose e = H'_/(/jfimaxg(2) |B(x)|), fc = 2W-. Then </>(x) > 1 and 

r°° fit r°° rli f°° tit 

/ e-(«"+^(*-0(a:))|<y e" At (i - 1)| < jf e"* = < | log + 0(1) 



as A — ► 0. Furthermore, on integration by parts 



'" 00 „ ks ds . e- fc 2 



implying that 

(79) ff (l)>e- fe /(fc + 2). 
It now follows from @, @, and @ that 

N(P W + A, Q)<(W + \)I-^r(k + 2)e k (\ log A" 1 ])^- 1 |Q( 2 > | max \B(x)\ 
~ 7m q(2) 

and (fr3|) is proved. □ 

The method of proving Lemma || also yields Proposition [| at the end of this 
section for an unbounded 0, and an operator P;y(f2) with a non-empty essential 
spectrum. However, we leave the result till then so as not to break the flow of the 
argument leading to the proof of Theorem ||. 

Lemma 5. Let Q = [0, R) 3 , W constant on Q, and \B\ > k > 0. Then, for any 
A>0, 

(80) lim{ A1 - 1 ^ 2 ^(P lv + A,g)--i 7 / \B\(W + X)ldx} =0 
h^o 2tH J q 

with 2fihn > W- . For any 7 > 

(81) lim{ A i^ 1 n 2 A/ 7 (P vv , Q)-/? 7 / \B\Wl + ^dx} = 
n -> Jq 

with 2fihn > in which 



Proof. As in [ pj| , Lemma 5.1], let f2^ , j G N, be non-overlapping squares obtained 
by translating = [0, R) 2 parallel to the coordinate axes to form a tesselation 
of M 2 . Choose a square C M 2 such that 

fih- 1 — [ B(x)dx = NeZ 

27T J D (2) 

with B extended by periodicity from Q (2 *> to R 2 . Let M = #{j : ttf } C £> (2) } and 
suppose, without loss of generality, that 

(82) (l-e)|^| < fi.fr 1 M\^\ < \N\, $ = — / B{x)dx 



2tt 



Q< 2 > 
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for e G (0, i). Set = x [Q,R) and denote by the torus with funda- 
mental domain £)(«*), d = 2,3. It is proved in [fH, (5.3)] that if [W + A]_ < 2fihn 

(83) limsup^ _1 /i 2 7V(Pvi/ + A, Q) < —~ [ \B\(W + \)ldx, VA > 0. 
h-^o 2?r^ J q 

For the lower bound we need A > in order to use Lemma ^. Let 

S(2)= xW\u*£ li ni(r»f ) ) 

S| 2) = {xGX( 3 ) : diai(a;, 5< 2 )) < 5} 

for 6 £ (0,f), and set 5^ = M 2 ) x [0,i?), S { s 3) = sf 3 x [0,i?). Then, JfW C 
u jii int(nj°) for d = 2,3. Also, since |B| > k 

5F*|S< 2) |< &S sm \B\dx 

= ^I D ^\ B \d^-^Mj ni2) \B\dx 



< 



e\N\ 



from (§3); note that £> is of one sign since B is continuous and \B\ > k. Therefore, 



and 



which implies that 



|S (2) I < -^re|AT| 
Kfin 

\SP\<IL^\N\ + M6 



|Sf J |< *-^\N\R + MRS 



(84) < M*\$\R + MR6 

— k{1— e) 1 1 

< C(e + <5)Mi? 

where C depends on n and $. We now proceed as in (|3^). The analogue is 

([F w+x (xW) + £f]fj) > (p w+x (sl 3) )M,M) 

for / in the domain of PwiX^) and a partition of unity {ipj}^ subordinate to the 
covering of by SP \jf =l nf } which satisfies (|l|) with K = M and 5 replacing 
pr. From this we conclude that 

(85) N (¥ w + \ + ^ 7 X^) < N(P W + A, Sjp) + MN(F W + A,Q) 
From Lemma || 

(86) N(F W + A + X&) > —R(W + A + - \N\ 

S z irn S z 

and by Lemma |I] and (Q) 

(87) N(¥ w + A, Sf } )<a^~ 2 | log A" 1 1 |Sf } |<^~ 2 | log A _1 |(e + 8)MR. 
It follows from (||)-© that 

1 (7^2 i 

MiV(P w + A, Q) > —R(W + A + - \N\ - C/i/r 2 | log A^Ke + i)M. 

Since 

\(W + A + - (W + X) 1 \<h/S 
o ~ 
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we have from ( |S2| ) 

H-WNfPw + X,Q)>-(W + A)£ |$| - C(-=.R|$| + + | logA^Ke + <5)i?) 

7T 



which yields (|80| ) on recalling (p3|) . Note that (|87| ) gives, for any 7 > 0, 

/>oo 

(88) M 7 (P W , Sf } ) = 7 / A 7 ~ 1 iV(P M / + A, 5j 8) )dA<At»~ a (e + 

Jo 

A similar argument then gives (pll). □ 



(89) lim ——{N(F W + A, Q) - ?T a «8 (pR|B|, W + A, Q)} = 



For \B\ > k and 2^hn > W- we have 2k(ih\B\ + W > 2^hn - W- > 0, for k > 1, 
and, for any 7 > 0, 

< B 1 ((j,h\B\,W,Q) = (ih^ [ \B\Wl + ^dx. 

Jq 

This fact gives 

Corollary 1. Let Q = [0,i?) 3 , VF a constant on Q, and \B\ > n > 0. Then, for 
any A > 0, 

I 3 

uniformly in fi > 0. Moreover, for any 7 > 

h 3 

lim — {MJ¥ W ,Q) - h- 3f BJuh\B\, W,Q)\ = 0. 

ft->0 fin + 1 

In order to remove the assumption that \B\ > k, we prove the following 
Lemma 6. Let Q = [0, i?) 3 and W a nonzero constant on Q. Then, for A > 
limsup^ A |iV(P w + A,g)-^ 3 » (^|B|,M/ + A,Q)| 
< CllogX-^WllQlmaxqmlBix)], 

and for any 7 > 

limBup^o-j&j |Af 7 (Pn,,Q)-n- 3 Q5 7 (^|B|,T^ + A,Q)| 
< CW2 +2 [|logVF_| + l]|Q|max Q(2) \B(x)\. 
Both ( |ff6[ ) and ZioZd uniformly in [i > 0. 
Proof. By Theorem [l] 

ft 3 

(92) lim sup — — -|JV(P W r + A,Q)-«8(/*fi|B|,W + A,Q)|=0. 

From Lemma |J 

-^-7V(P W + A, Q)<| log A" 1 !^ + A)i|Q| max |B(a;)| 

/i/l +1 ~ Q( 2 ) 

and 

jjg r n- 3 » (/ifi|S|,W + A,Q) 

< f Q \B\[(W + X)i + 2J2 k >i( 2 klJ i n\B\ + W + X)i]dx 

< Jq( W + x )h\B\ + {nK)- l {W + A)_]dx. 
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Therefore, 

su P/jfi > c j^WiPw + A, Q) - h- 3 <B(fih\B\, W + X,Q)\ ^ 
<|logA- 1 |(T^ + A)^|Q|max Q ( 2) \B(x)\ + ^\Q\{W + A)! 

and (^0|) follows since C is arbitrary. The proof of ( |9l| ) follows in a similar manner. 

□ 



Theorem 6. Let B(x) = (0,0, B{x)). Suppose that 

1. B is continuous, 

2. W- € L°°(Cl), 

3. for some p > l,W^_,bpW- S L 1 (f2), where b p is defined in ([7 
Then, for all 7 > and A > 0, 

(93) lim — \M^(F W + \,n) - h- 3< B-«Uh\B\,W + \,n)\ =0. 

ft->0 fin + 1 1 J 



uniformly for p>0. If 7 > 1/2, ZioWs /or A = 0. 

Proof. We prove the result for 7 = 0, the proof for 7 > being similar. The 
last assertion in the theorem, concerning 7 > 1/2, is proved i n p3| . Let Wo, Qk, 
k = 1,2,..., k be as in §2 (Ai)-(A 3 ). As in §6 of Sobolev |f7for 77 > 0, we 
partition each Qk m to a finite number of cubes Qjk, j = 1,2,..., such that for 
each pair k, j 



(94) \B(x)-B(y)\<n, x,yeQ k 



3 ' 



Let J+ = I+{rf) := {(k,j) ■ rnaxQ fc . |-B(x)| > 2r]} and 7_ = the complemen- 

tary set. In view of ( p^ ) we have that |-B(a;)| > r\ for x 6 Qfcj, (k,j) € i+. Since 



A' 



7V(P Wo + A, SI) > £ iV(P Wo + A, Q k ) 
fc=i 

and from Corollary [l] and Lemma ^| 

liminf ri ^ ^ {iV(P Wo + A, 0) - /r 3 <B (^|B|, W + A, fi)} 
> -^E (Me /_ / Qfc , 3 (W + A)*ds> - 77. 

Since 77 is arbitrary, then 
h 3 

(95) limmf— — - \N(F W + X,fl) — h~ 3f B (ph\B\, W + A, SI)} >0 
R^o /i/l + 1 

for each A > 0. For the reverse inequality we proceed as in the proof of (p7|). Let 
the interior of each Qk be denoted by int(Qk). Set 

. , S := Cl\u£ =1 int{Qk), S pr := {x e SI : dist(x, S) < pr}, 

^ ' Q{pr) := {x e SI : dist(x, Q) < pr} 
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for some p £ (0, 1). Construct a partition of unity {ipk}k=-i subordinate to the 
covering Li k=1 int(Qk) U S pr of £1: 

(97) 

(*) fier(Sl), ipo&C$°(S pr nQ pr ), ip k £C%°(int(Q k )), k = l,...,K; 

(H) Efc=-i^f = l.Va;; 

(m) V-iC 1 ) = 0j x e Q, tp-i(x) = 1, x £ Q(pr); 

(*«) Ef=-i |V^fe(2;)| 2 < C(pr)- 2 x Spr nQ( P r) 

for some C > where Xpr is the characteristic function for S pr D Q(pr). It follows 
that 

(98) |5 pr nQ|<p|Q|, |Q(pr)\Q|<pr, 

and for every / e [C£° (ft)] 2 

£ PaWM) < (VofJ) + Ch 2 (pr)- 2 ( X s pr n Q ( P r)fJ)- 
fc=-i 

Let r = Aft and A - C^p)" 2 > 0. Then, 

m (CV + A)/, /) > Ef=-i(PPwb + A - Cft 2 (pr)- 2 Xpr ]V fc /, Vfc/) 

> Ef=o([ F ^ +A-C*ft 2 (H" 2 Xpr]^/>fc/) 

since V'-i = in Q,and Wo = for x $ Q. This implies that 
n s N(P Wo +X,n)< N(P Wa +X-Ch 2 (pr)- 2 ,S pr nQ(pr)) 

[ ' + Ef=i # (HV + A - CT^r-r 2 , Qfc)- 

From Corollary |l| and Lemma || 

limsu P^o7ifeEf=i {N(P Wo +\-{^,Q k ) 
(101) -ft- 3 » (^|B|,Wo + A-^,Q fe )} 

< ^Ecm^/q^C^o + a)!^. 

and from (§|), 

nn9 , ^|®o(^|B|, W + A - gr, Q k ) - ® o 0*>i|B|, W„ + A, Q k )\ 

^Ap f (A7F- 



From ([73|) and fl9J) 

Ai - 1 ft 2 7V(P H / + A - CH 2 (pr)- 2 , 5pr n Q{pr)) 
1 l 



(103) <eV(^) 2 i p |Q|. 



It now follows from ( IOC )- (|103| ) , on allowing h — > 0, A — * oo, 77 — > 0, p — > 0, in that 
order, that 

ft 3 

limsup— — -{N{¥ Wo + A,fi) -ft- 3 «B (^|B|,VK + A,fi)} < 
n^o /xfi + 1 

which together with ( |95| ) implies 
ft 3 

(104) lim — — -{N(¥ Wo + A,n) - ft- 3 » (^|B|, W Q + A,0)} = 0. 
i-o fin + 1 
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Similarly, we have for any 7 > and A > 

h 3 

(105) lim — — -{MJV Wo + A, 0) - /T 3( 8 7 (^|B|, W + A, VI)} = 0. 

h^o /in + 1 

The proof for general W is similar to the corresponding part of the proof of The- 
orem ^| in §3, using the Weyl inequalities in the same way as Sobolev does in p3| , 
but, instead of ( E6I) w e now use the inequality derived by Shen in Lemmas 5.1, 5.3, 
5.4 and (5.14) of H20|, namely, 



(106) H 3 N(¥ W + \,R 3 )<^{ / \W-(x)\ 2 dx + iih / 6 p (x)W_(x)dx}. 

~ V A JR3 J R 3 

□ 

The following analogue of Theorem |^ is also readily proved. 

Theorem 7. Let B(x) = (0, 0, B{x)). Suppose that 

(i) W , B are continuous, 

(ii) |W(x)| — * uniformly as |x| -^00 in fi. 
Then, @ holds for all 7 G [0, 1) and A > 0. 

Finally in this section we give the result promised after Lemma 

Proposition 5. Let fl = fl^ x R, where is of finite measure in R 2 , and 
suppose that 

(i) B = (0,0, B(x)) and |B| S L°°(n^); 

(ii) infn W(~x) > —cx, R R) (^3), x = (x, x 3 ), for some positive numbers c, R. 
Then for A > 

(107) 7V(P H / + A,f7)<^ 3 (^+l) v ^|logA" 1 |i? / (\B{x)\+c)dx 

~ JO< 2 ) 

and, for any 7 > 0, 

(108) M 7 (Pvi/,^)<^ 3 (^+l)c 7+ ^(|logc| + l)i? / {\B{x)\+c)dx. 

Proof. In the notation of the proof of Lemma |], 

Pw > Po 2) + M 2 <9| - W) 

where W(x) = W(x 3 ) = c\ { _ R R) (x 3 ). The expression —h 2 d 2 — W has a self-adjoint 
realization in L 2 (R) with essential spectrum [0, 00) (see jl^], Proposition 2.1), and it 
is readily shown that there are negative eigenvalues at the solutions of the equation 



(109) tanffv^A)=4 ^ A( " A) , Ae(0,c) 

n c — Za 

corresponding to the eigenfunctions 

(110) 

(p{x 3 ,X) = 

e-$*», x 3 >R, 
e^ Xs , x 3 < -R, 
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where 6 := — 1. If these negative eigenvalues are denoted by — e m , m 
M, then e m < c and 



(in) 



M<#{n:c-(^!>0} = 



We now proceed in a similar way to the proof of Lemma [|. We have that for any 
£€(0,1) 



(112) 



N(P W + \,n)< J2^ =1 N(P { 2) + A - e mi nW) 

< 2 E™=i N(eH^ 2) - efih\B\ + A - e m , Sl^). 

Now substitute 

^l = e^ 2) +A, B = -e7i 2 A + A, V = {e»h\B\ + e m )x n(2) 

in Theorem 2.4 of [fill. This gives, as in @, 



N(A— V) < 



V{x)dx 



,-At 



4>(x) 



(t -</>(*)) 



dt 



where <j){x) = k/V(x) and g(l) satisfies (78|). Choose 

e= , A; = max + c. 

1 

Then 0(a;) > 1 and the result follows as before. 



□ 



6. The Dirac operator: magnetic fields with constant direction. 
Theorem 8. Suppose that 

1. B is continuous, 

2. V- e L°°(fJ), 

3. /or some p>l, V*, b p V* € L^O), 

4. |{x€fi: |V(x)| > 2}| = 0. 

Then, for all 7 > 0, 

Um A ^ + limft^o ( 7 £ T ){i[M y (B y , fi, 1(A)) + Af 7 (D_ v ,n,/(A))] 

-h- 3 <B 7 ( t ih\B\,\ 1 ,n)} = 0, 

w/iere Ai = 2\V\ + V 2 and I(X) = (-V1 - A, \A~~X). 

Proof. The proof follows from Theorem]^ just as Theorem]^ follows from Theorem |[ 
As for (55) and (|60|), we show that (in the notation of (Ai)-(A^) in §2) 

liminf r ^ (_*L.){(l/ 2 ) [JV((B Vo) 0, /(A)) + 7V((D_ yo , O, /(A))] 

-h- 3 [*B(jJi\B\, -A?, n) + <B(nH\B\, -A°!,n)]} > 



where A%(A) = (v 7 ! - A ± |F |) 2 - 1, and 

limsup^o ( ? g T ){(i/2) [JV((D Vo , n, /(A)) + AT((B_ Vo , O, /(A))] 

-^- 3 [?B(/xn|B|,-A?,0)-Q5(Mn|S|,-A° 1 ,0)]} < 

The rest of the proof is similar to that of Theorem ||. □ 
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7. Appendix 

7.1. Approximation of gauges in cubes. Let XQ fc = (xg k , Xg k , Xg k ) be the 

center of the cube Qk- Since div(B) = 0, we may choose the gauge a — (a%, a 2 , 0) 

ai(x) := B 3 (xi,t,x^ )dt + f*f B 2 (xi,x 2 ,t)dt, 

a 2 (x) := i/^ 1 B 3 (t,x 2 ,x%)dt- B 1 (x 1 ,x 2 ,t)dt 

Qk Qk 

- see the proof of Lemma 3.1 of Sobolev |25|] . Let B°(x) := B(xQ fc ), x£ft and 
define the piecewise linear vector potential a by 

ai(x) := -^B 3 (x Qk )(x 2 ~x 2 Qk ) + B 2 {x Qk )(x 3 -x 3 Qk ), 

1 , , 

a-2(x) := -^(xqJ^i-XqJ-B^xqJ^s-XqJ, 
a 3 := 0. 

for x e Qfe. Then V x a = B° in Qk and (11) is satisfied. 
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